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Within the s-d model description, we derive the current-driven spin torque in a ferromagnet, 
taking explicitly into account a spin-relaxing Caldeira-Leggett bath coupling to the s-electrons. We 
derive Bloch-Redfield equations of motion for the s-electron spin dynamics, and formulate a gradient 
expansion scheme to obtain non-adiabatic (higher-order) corrections to the well-known adiabatic 
(first-order) spin torque. We provide simple analytical expressions for the second-order spin torque. 
The theory is applied to current-driven domain wall motion. Second-order contributions imply a 
deformation of a transverse tail-to-tail domain wall. The wall center still moves with a constant 
velocity that now depends on the spin-polarized current in a non-trivial manner. 
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I. INTRODUCTION 

Many recent experiments have brought to attention 
the intricate and rich interplay between spin-dependent 
transport and the magnetization dynamics in ferromag- 
netic (FM) materials [l|, [H, Q . A particularly interesting 
issue concerns the spin torque T 0, @, @, H, ex- 
erted by conduction electrons on the FM magnetization, 
—M s n, where n(x, t) is a space-time dependent unit vec- 
tor and M s is assumed constant. The spin torque is cru- 
cial for the understanding of many topics of present scien- 
tific and technological interest, including current-driven 
domain wall (DW) motion [lOj, magnetization preces- 
sion and switching in multilayer geometries [TT| , and spin 
transport in general [l2j. Here, we study the spin torque 
exerted by a steady-state dc spin current J s [for the def- 
inition of J s , see Eq. (|10[) ] in the presence of spin relax- 
ation. Common microscopic mechanisms responsible for 
spin relaxation are, for instance, provided by magnetic 
impurities or the various sources for spin-orbit scatter- 
ing, but we shall focus on generic features and study 
a phenomenological Caldcira-Leggett-typc heat bath de- 
scribing random magnetic fields with Gaussian statistics 
which are responsible for spin relaxation. 

As has been discussed recently, spin relaxation is es- 
sential in determining the non-adiabatic corrections to 
the well-known adiabatic spin torque [H, Q 

T ad = -J s d x n. (1) 

It is common practice to call all terms beyond Eq. ((T|) 
'non-adiabatic' (but see below). In particular, Zhang and 
Li [l3| predicted a leading non-adiabatic spin torque con- 
tribution 

Tzl = PJ s n x d x n, (2) 

which was confirmed by later work [3, EH, fl6|. [TtI [l8| , 
although some questions were raised in Ref. Here 
the dimensionless (3 parameter is given as (3 = fi/Ar s , 
where A is the FM exchange splitting and t s denotes a 
spin relaxation time. Although (3 is generally small, Tzl 



can profoundly influence many of the above-mentioned 
phenomena. Recently, a complicated and highly non- 
local form of the spin transfer torque including non- 
adiabatic corrections due to fast varying magnetic tex- 
tures has been proposed [2(| ■ This form includes indirect 
Ruderman-Kittel exchange interaction processes 
between the localized spins mediated by itinerant elec- 
trons. A scattering approach has been employed in Ref. 
[23l | to calculate the spin torque in the limit of a sharp 
DW. However, it is difficult to deal with spin relaxation in 
that formalism. Strong enhancement of the spin torque 
and of the domain wall mobility due to spin-orbit cou- 
pling in magnetic II- V semiconductors has been reported 
in Ref. Q. 

In this work, we provide a gen eral derivation of the spin 
torque within the s-d model po1.[26|. see Sec.|TTJ The self- 
consistent Stoncr model gives very similar results for this 
problem but is technically somewhat more demanding 
HUGH E3 . While spin torques in the presence of spin re- 
laxation have so far been studied mostly for small pertur- 
bations around a homogeneously magnetized ferromag- 
net [HI, EE EH , we here describe the theory for the lead- 
ing nonadiabatic corrections for non-homogenous (both 
in space and time) magnetization profiles. For these cor- 
rections and not extremely large applied currents, the 
s-d model is expected to provide a reasonable descrip- 
tion. We note in passing that recently a simple classi- 
cal model of non-adiabatic current-induced spin torques 
has been proposed [23] ■ In our work, spin relaxation 
is included within the weak-coupling Bloch-Redfield ap- 
proach [H, [29|, HH, where the s-electrons are cou- 
pled to a spin-relaxing environment modelled by a phe- 
nomenological harmonic oscillator bath (32j, see Sec. IIIII 
Moreover, the itinerant electron spins are influenced by 
a space-time dependent magnetic field produced by the 
localized electron spins. This external field also deter- 
mines the stationary state to which the itinerant spins 
relax. Validity of the Bloch-Redfield approach requires 
weak coupling of the s-electrons to the bath, i.e., suffi- 
ciently long spin relaxation times t s , and a bath memory 
time t c short compared to all other relevant timcscalcs. 
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Both conditions hold in most cases of practical interest. 

Based on this formalism, we then determine the spin 
torque T acting on the FM magnetization n(x, t) within 
a gradient expansion around the adiabatic limit, see 
Sec.|lVl We illustrate the method by computing all first- 
and second-order derivative contributions to the spin 
torque, i.e., all terms that involve d X}t n, d x t n, d x d t n, or 
products of first-order derivatives. The first-order terms 
provide the leading terms in this expansion, which we 
call the adiabatic torque. Note that our convention dis- 
agrees with common usage but emphasizes that these 
terms completely dominate the spin torque for slowly 
varying n(x,t). Indeed, a first-order expansion recovers 
both Eq. (p} and |J2J). This shows that the 'non-adiabatic' 
term ^) is really a consequence of spin relaxation rather 
than being related to higher-order terms in the gradi- 
ent expansion [3, EH, EB[. Our approach then allows 
to systematically evaluate higher orders in an expansion 
around the adiabatic limit. It should be stressed, how- 
ever, that the s-d model may become difficult to justify 
for extremely sharp or fast features in n(x,t) (e.g. very 
narrow domain walls), and in practice this expansion is 
only useful for the first few orders. 

Here we will focus in detail on the second-order spin 
torque contributions which turn out to be of relatively 
simple form. To illustrate some physical consequences of 
these new terms, in Sec.[V]we address the problem of DW 
motion driven by a spin-polarized current for the case 
of a tail-to-tail transverse DW. We find that the second- 
order contribution leads to a deformation of the adiabatic 
DW shape. Nevertheless, under steady state conditions 
the DW center moves uniformly with constant terminal 
velocity Vdw- Importantly, this velocity depends non- 
trivially on the spin current density J s . After following 
the known linear behavior of the velocity found from the 
first-order calculation for small J s , Vdw does not grow 
further with J s but instead starts to decrease. We expect 
that when including all orders, this would indicate sat- 
uration behavior. Our results seem in qualitative agree- 
ment with experimental observations, where measured 
DW velocities are typically smaller than those predicted 
by first-order calculations (33[. Moreover, the deforma- 
tion of the DW profile is also observed in experiments and 
can even induce a transformation to a completely differ- 
ent DW type [Hj]. While the quantitative description 
of experiments on current-driven DW motion is likely 
to require more refined models, our principal aim is to 
provide general expressions for the leading corrections to 
the known spin torque expressions dTJ) and ([2]), and to 
illustrate typical effects caused by these new terms. The 
practical importance of such effects was also revealed by 
recent numerical studies [HI, l3rl [37j . 

Our paper then closes by offering some concluding re- 
marks in Sec. IVIl The derivation of the relaxation kernel 
under the Bloch-Rcdfield approach can be found in the 
Appendix. Throughout the paper, we put h = ks = 1. 



II. MODEL 

We consider an (infinitely long) FM wire with homoge- 
neous magnetization along the cross-section, n = n(x, t), 
where x is the longitudinal direction. This simplifica- 
tion allows us to work with an effectively one-dimensional 
(ID) theory. For thin FM nanowircs with only a few oc- 
cupied transverse sub-bands, this description is directly 
relevant but for wider wires, nontrivial transversal 
magnetization profiles (such as vortex walls) are thereby 
excluded. Moreover, the ID model also neglects spin 
waves propagating in the transverse direction that can 
be excited in wider wires. We note in passing that in 
strictly ID metallic wires, electron-electron interactions 
can cause non-Fermi liquid effects [39[ that could also in- 
fluence domain-wall motion [5(|. However, such effects 
are not studied in what follows. 

Adopting the widely used s-d model for itinerant FMs 
@, B HE , the relevant dynamical degrees of freedom 
are (i) the FM magnetization n(x, t) associated with lo- 
calized d-electrons, and (ii) the conduction electron spin 
current density 3(x,t) and spin density s(x, t), respec- 
tively, describing the dclocalized s-clcctrons. The FM 
magnetization then obeys the Landau-Lifshitz-Gilbert 
(LLG) equation Q 

d t n = -7 n x B off [n] + an x d t n + T, (3) 

where 7oB e fj includes external, anisotropy, and ex- 
change magnetic fields unrelated to the coupling to s- 
electrons. The Gilbert damping parameter a provides a 
phenomenological description of dissipative influences on 
n (again in the absence of s-electrons) . Our principal 
aim is to find the spin torque T when the FM carries 
a (spin-polarized) current and both spin relaxation and 
magnetic texture (spatio-temporal variations of n) are 
present. 

Let us then address the conduction electron degrees of 
freedom. Expressed in terms of right- and left-moving 
(p = R/L) quasi-particle annihilation operators c PjCr (x) 
for spin a =1,1, the ID spin and spin current density are 
defined as 

s{x)=J R (x)+3 L (x), J(x)=v(3 r (x)-J l (x)), (4) 
where v is the spin velocity and the chiral currents are 

3r/l{x) = - : c R/L {x)ac R/L (x) : . (5) 

Here colons denote normal ordering, spin indices are left 
implicit, and a are standard Pauli matrices acting in spin 
space. The currents ([5]) obey the SU(2) Kac-Moody al- 
gebra [3i|. Within this ID description, the low-energy 
Hamiltonian describing the conduction electrons is uni- 
versal and given by 

H = -iv ^2 P dx c\d x c p . (6) 

p=R/L=± •* 
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The spin sector of Hq can now be expressed directly in 
terms of the currents ([5]), and completely decouples from 
the charge sector. (In a strictly ID wire, this procedure 
works even when forward-scattering electron-electron in- 
teractions are included (39[.) The resulting Sugawara 
representation of the spin part of Hq (the charge part is 
irrelevant and omitted here) is 



H = 



dx : s • s 



1 



r J J 



(7) 



Within the s-d model, s-electrons are coupled to the 
time-dependent magnetization n(x, t) only by exchange 
processes. Effectively, the magnetization then acts on the 
s-electrons as a time-dependent external field, 



H ex (t) — A / dx s(x) • n(x,t), 



(8) 



where A is proportional to the FM exchange coupling 
and can have either sign. For the resulting 'system' (the 
s-electrons), we then arrive at the Hamiltonian Hs(t) = 
Hq + H ex (t) . Equation ([8]) now generates a spin torque 
entering Eq. ([3]), 



T = n x - 



SH ex 



= —An x s. 



(9) 



Our task is then to compute the spin density s(x, t) in the 
presence of an external dc spin current J s and for a given 
dynamical magnetization profile n(x,t). We mention in 
passing that we do not perform a unitary transformation 
into an adiabatic reference frame where the local quan- 
tization axis for the itinerant electrons is aligned with 
n(x, t). Albeit this transformation has been used in many 
papers on the subject, see e.g. Ref. and references 
therein, we believe that it docs not offer advantages in 
the presence of spin relaxation. 

In order to enforce a finite spin current, in principle 
(i) we need to resort to a non-equilibrium formalism. In 
addition, (ii) we need a proper description of the relax- 
ation processes driving the momentary state to a station- 
ary state. Let us first discuss issue (i), where a simple 
boundary condition can be used instead of the full dy- 
namical formalism, see also Ref. (38J. We want to de- 
scribe a steady-state situation with externally imposed 
constant spin current 



J =™ 



(10) 



expressed in terms of the charge current I flowing 
through a wire of cross-section A with spin polariza- 
tion < P < 1. Since this stationary spin current is 
externally enforced, we include it as a boundary con- 
dition in the equation of motion of the itinerant elec- 
trons derived below. This boundary condition has to 
be imposed far away from any magnetic texture (i.e., 
at x — > Too) j where it implies a constant spin current 
J = J s n. Only around, say, a DW center, the true 



spin current density deviates from J(x,t) = J s ii(x,t). 
Coming to issue (ii), it is then tempting to identify the 
steady state to which an actual spin current configuration 
tends to relax with J(x, t) = J s xi(x,t). This indeed we 
find, see Sec. Mil The Heisenberg equations of motion, 
dt^n/L = i[Hs(t), «Jr/z,]j now yield operator equations 
for the spin density and the spin current, 



d t s + d x J 
d t 3 + v 2 x s 



-As x n + T s , 
-AJ x n + r, 7 . 



(11) 
(12) 



The relaxation terms r s /j[s, J] due to the coupling of the 
s-clectrons to the spin-relaxing environment are specified 
in Sec. IIIIl Far away from magnetic textures, both the 
derivative terms and the relaxation terms T Sj j are irrel- 
evant, and taking s = son and J = J$n with arbitrary 
constant coefficients so and Jo solves these equations. 
The boundary condition on J discussed above then en- 
forces Jo = Js- To determine s , one has to find the 
stationary state of the equation of motion (jlip in the 
presence of spin relaxation. We here anticipate the re- 
sult so = — X S A derived in Sec. IIIIL see Eq. ([TBI) , where 
Xs = (2ttv)~ 1 is the spin susceptibility. 

To summarize our discussion, we seek solutions to 
Eqs. dTTJl and {T2J) of the form 



3{x,t) = J s n(x,i) + ^J fe (x,t), (13) 

fe=i 

s(x,t) = -( Xs A)n(x,t)+J2 s k{x,t), (14) 

k=l 

where the fcth-order terms contain fcth-order spatio- 
temporal derivatives of n, or respective products of lower- 
order derivatives. Such an expansion provides a system- 
atic way to classify deviations from the adiabatic limit 
where only the k = 1 terms are retained. However, as we 
have discussed above, in practice only the first few orders 
in this expansion can be reliably extracted under the s-d 
model description. 



III. SPIN RELAXATION 

For a description of the spin-relaxing environment act- 
ing on the spin dynamics of the s-electrons, we employ 
a Caldeira-Leggett system-bath approach [32| , where the 
s-electron spin density and spin current are linearly cou- 
pled to fluctuating magnetic fields B Sj j(x), 

H SB = j dx [B s (x, t) ■ s(x) + Bj(x, t) ■ J(x)]. (15) 

These fields provide a simple phenomenological mod- 
elling of the effects of random magnetic impurities and of 
spin-orbit scattering. The Ansatz (p~5|) assumes 'soft' for- 
ward scattering due to the fluctuating fields, where right- 
and left-moving fermions retain their chirality. Back- 
ward scattering processes, for instance associated with 
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clastic potential scattering [38|, will modify our quan- 
titative conclusions, as discussed in more detail below. 
Note also that the fluctuating fields are assumed to 
solely couple to s-electrons, and the FM magnetization 
n(x, t) is only affected indirectly via its exchange cou- 
pling to the s-electron spin density s(x,t). The Hamil- 
tonian Hg describing the uncoupled fields B s .j corre- 
sponds to a suitable bath of harmonic oscillators [32| . 
With (B s .j) = 0, where (O) denotes the equilibrium 
quantum-statistical average of some operator O with re- 
spect to Hb, all relevant properties are thus encoded by 
specifying their two-point correlation functions. The to- 
tal Hamiltonian is then H to t(t) = Hs(t) + Hsb + Hb, 
with H s (t) = H + H ex (t) as given in Eqs. |7J) and ©. 

To describe the effects of Eq. (|T5|) on the s-electron 
spin dynamics, we adapt the textbook weak-coupling 
Bloch-Redfield approach [H, [U, HO, [Ml to a field theory 
as required here. For the convenience of the interested 
reader, we provide a detailed discussion of the Bloch- 
Redfield equations in the Appendix. The basic assump- 
tions underlying this approach are (i) short bath corre- 
lation times, i.e., the two-point correlation functions of 
B s .j(x, t) vanish on a timcscalc r c small compared to all 
relevant other timescales, and (ii) weak system-bath cou- 
pling Hsb- Note that Hsb in Eq. (fl"5|) is of the general 
form considered in the Appendix, see Eq. (| Al|) . This im- 
plicitly also requires that we stay not too far away from 
the adiabatic limit. Following the steps in the Appendix, 
we obtain the relaxation term entering Eq. in the 
form 



a(x,t) + (x s A)n(a:,t) 



(16) 



where \s — (2irv) and t s is the spin relaxation time. 
Here r s is assumed uniform in the x, y, z directions; the 
generalization is discussed in the Appendix. Since n(x, t) 
acts effectively as a space-time dependent magnetic field 
on the spin density s, see Eq. it directly appears in 
the stationary state approached in the relaxation term 
(|16[) . In a similar manner, we obtain 



J(x,t) — J s n(x,t) 



Tj 



(17) 



For the simple system-bath coupling (fT5|) . which disre- 
gards backscattering processes, the relaxation times r s 
and tj turn out to be equal. In the following, we shall 
allow for different relaxation times tj ^ t s , which re- 
sults when allowing for more general models involving 
backscattering processes [HI, i.e., when elastic disorder 
is present. 

Before discussing the resulting spin torque, let us 
briefly comment on the validity regime of our descrip- 
tion. Both the derivation of the relaxation terms, see 
Eqs. (fTB|) and (jTTJ) , and the iterative procedure for solv- 
ing the resulting equations, see Sec. IIV1 suppose that 
we are not too far away from the adiabatic limit. By 
that we mean that the space (time) variation of the mag- 
netization direction n(x,t) is sufficiently smooth (slow) 



when compared to a characteristic lengthscale A (time 
scale X/v). This 'spin transport lengthscale' has been es- 
timated by Zhang and Li, see Ref. [l3|]. Quantitatively, 
we suppose that both dimensionless non-negative func- 
tions X\d x n(x,t)\ < 1 and (X/v)\d t n(x, t)\ <C 1. Un- 
der these conditions, an iterative expression like Eq. Q18p 
implies a well-behaved perturbation scries, where higher 
and higher orders give smaller and smaller corrections. 



IV. SPIN TORQUE 

We next discuss how to extract the spin torque from 
an iterative solution of Eqs. (jTT|) and (fT2"|) . To that end, 
we perform a gradient expansion indicated by the series 
in Eqs. (fT3"|) and p4[) . This gradient expansion then gen- 
erates a series for the spin torque, 



T(x,t) 



k=l 



(18) 



Of course, the zero-order terms in Eqs. (jT3|) and Q14p do 
not generate a torque. This expansion naturally suggests 
to define the adiabatic spin torque as Ti, the dominant 
piece for sufficiently slow and smooth variation of n(x, t). 
As remarked above, this definition differs from standard 
usage but provides a systematic way to classify different 
spin torque contributions. 



A. Adiabatic spin torque 

Let us now compute the first-order spin torque Ti. 
The derivatives of Si or Ji are of second order and can 
be dropped here. Inserting the Ansatz (fT3"|) and (fT4")) into 
Eqs. (IT]) and (fTJ]) then gives 

(t^ 1 -Anx)si = x s A<9 t n - J s d x n = a x , 
(r7 1 -Anx)J 1 = -J a d t n + v 2 XsAd x n = b 1 .(19) 

Here Ji and Si completely decouple. Exploiting that the 
equation Ax + y x x = v with real A and real 3-vectors 
x and y is solved by 



A 2 v + Av X y + (v ■ y)y 
A(A 2 + |y| 2 ) 



it is straightforward to obtain 



Asi = -ai x n + /5ai, 
AJi = -bi x n + Phi, 



(20) 



(21) 
(22) 



where we have introduced the dimensionless parameters 

(3 = H/(At s ), P = h/(A TJ ). (23) 

In accordance with the weak-coupling Bloch-Rcdficld ap- 
proach, in numerical prcfactors we have neglected terms 
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of order (3 2 and (3 2 against unity. The parameters j3 and 
(3 are small expansion parameters in this approach. 

For the discussion of the adiabatic spin torque Ti, we 
only need si as given in Eq. (|21[) . Inserting this into 
Eq. ([9]) gives the complete first-order spin torque, 

Ti = - J s d x ii + /3J s n x d x n + A\ s [d t n - a'n x <9 t n] , 

(24) 

where a' = [3. The first term is the adiabatic spin 
torque JT]), and the second term is the 'non-adiabatic' 
correction ([2]). Moreover, the third term can be ab- 
sorbed by renormalizing the gyromagnetic ratio 70 and 
related parameters, cf. Ref. [13(. Finally, the last term 
gives an additional contribution to the Gilbert damp- 
ing that is absorbed by a redefinition of a. Our as- 
sumption of equal longitudinal and transverse relaxation 
times, T XtV = t z = t s , cf. Appendix, strictly implies that 
a 1 = f3. Under more general conditions, however, this 
relation will not hold anymore The approach 

presented here can straightforwardly be generalized in 
this direction. Equation (|24[) thus confirms and repro- 
duces the results of Zhang and Li [l3j], and effectively 
yields the first-order (adiabatic) spin torque in the form 

Ti = Tad + Tzl- (25) 

B. Non-adiabatic spin torque 

Next we calculate the second-order terms in Eq. (fT5T) . 
The second-order contributions S2 and J 2 are obtained 
by inserting the Ansatz ([13]) and HU) into Eqs. (fTT]) 
and (Ti"2"|) , thereby exploiting Eq. (TT9"]) . Dropping terms 
of higher than second order, the resulting equations can 
be solved again, with the result 

As 2 = a 2 x n /3a 2 , (26) 
AJ 2 = b 2 x n-^b 2 , (27) 

with the auxiliary vectors 

a 2 = ftBi + flxJi, (28) 
b 2 = d t 3x + d xSl , (29) 

where again terms of order j3 2 were discarded against 
unity in prefactors. We have also omitted terms oc n 
resulting from Eq. (|2H)) since these do not generate a spin 
torque. In fact, only s 2 in Eq. (|26p is required. After 
straightforward algebra, we get 

T 2 = n x (c tx d t d x n + c tt d 2 n + c xx d 2 n (30) 

+ n x [d tx d t d x n + d tt d 2 n + d xx <9^n]J 

with the coefficients 

c tx = 2J S /A, c tt = c xx /v 2 = -Xs, (31) 

d tx = -{30 + 0)J s /A, 

d tt = 20Xa, d xx = v 2 X s([3 + 0). 



We stress that as long as the magnitude of the magneti- 
zation vector is conserved, n 2 (x,i) = 1, Eq. (f3T))) repre- 
sents the most general second-order spin torque allowed 
by symmetry constraints. At this stage, let us compare to 
the results of Ref. (l5| , where a spin torque contribution 
T' 2 = (J s /A)n x dtd x n was reported. This should match 
our c^-term, but differs by a factor 2. The reason for 
this difference is explained by noting that Ref. [l5[ em- 
ployed the self-consistent Stoner description, where one 
finds indeed a factor 2 difference to our s-d model pre- 
diction [13] ■ The other contributions to the second-order 
torque in Eq. (|30p . in particular the d te -term, have not 
been reported before. Interestingly, the spin torque con- 
tains J s -indcpendent and linear-in-J s terms only. How- 
ever, even the linear-in-J s spin torque (|30p leads to a 
nonlinear response of the DW velocity as a function of 
J s , see Sec. [V] 

The second order torque T 2 in Eq. (|3H)l introduces sev- 
eral new features into the LLG equation ([3]). First, the 
only terms proportional to the spin current J s are due 
to Ct x and d tx . All other terms are less important when 
large spin currents are applied. For instance, the c xx - 
term can be combined with the exchange term in the 
effective magnetic field B e g entering the LLG equation 
©, see below in Eq. (f3"3"|) , i.e., it can be absorbed by a 
renormalization of the exchange constant J. Moreover, 
the Cu and <itt-terms introduce acceleration terms in the 
LLG equation. Since we are interested in current-driven 
spin torques here, we will exclusively focus on the Ct x 
and dtx terms from now on. We have checked explicitly 
that the results below do not qualitatively change if the 
d xx , Cu, dtt terms are included in a perturbative scheme. 
The second-order spin torque is thus taken as 

rp 2J S / 3/3 +1 \ 

T 2 = — n x d t d x n — n x d t d x n I . (32) 

This is a central result of our work. Albeit the coeffi- 
cient dt x is smaller by a factor « 2(3 compared to Ct x , it 
can qualitatively influence the magnetization dynamics, 
just as is the case for the first-order terms. We believe 
that Eq. (|32p should be included in micromagnetic simu- 
lations of (for instance) DW dynamics. Below we will see 
that these new terms lead, among other phenomena, to 
a deformation of the shape of a DW, as the spin torque 
now involves the time derivative of the gradient of n. 

V. DOMAIN WALL MOTION 

In order to illustrate the effects of the new spin torque 
terms (|32[) . we now consider the problem of current- 
induced domain wall motion. Let us take a tail-to-tail 
DW, being a solution of the stationary LLG equation © 
with 

7oB cff [n] = Jd 2 n + Kn z e z - K ± n y e y , (33) 

where J is the ferromagnetic exchange coupling, K the 
longitudinal and K± the transverse anisotropy constant 
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(incorporating the effects of demagnetizing fields). The 
resulting initial DW configuration at time t = (when 
the spin current is switched on) with arbitrary center Xq 
is then given by 



n(x, 0) = cosh 



x — Xq 



w 



tanh 



x 



(34) 

where w = \f~JJ~K is the initial width of the wall. We 
then determine the time evolution of this configuration 
from the LLG equation in the presence of a spin polarized 
current J s switched on at t = 0, taking the full spin 
torque T = Ti + T 2 , see Eqs. (25) and (32). 

Due to the complexity of the resulting LLG equation, 
we have to solve it numerically. To that end, we spa- 
tially discretize the LLG equation into a discrete ID 
spin chain, with spins at x = ia (integer i) for lat- 
tice constant a. The resulting first-order ordinary dif- 
ferential equation (in time) is then solved numerically 
starting from the initial configuration (33). To be spe- 
cific, we choose a spatial grid with x £ [0, f 000a] and 
scale time in units of to = a/v (and energies corre- 
spondingly). In these units, we take [33] J = 3, K = 
0.1, K x = 0.01, A = 0.2, and x = 500a. This implies 
the initial DW width w ~ 5.5a. Moreover, we choose 
a = j3 — (3 = 0.03. Let us briefly estimate the corre- 
sponding dimensionful parameters. For this, we take the 
Fermi velocity of permalloy, v — 2.2 x 10 5 m/s [43|, with 
a unit cell of volume Vb = a 3 = (1 nm) 3 . The time scale 
to = 4.5 x 10~ 15 s follows. The values for the magnetic 
constants are J = 4.4 x 10~ 10 J/m, K = 14.5 x 10 6 J/m 3 , 
K x = 14.5 x 10 5 J/m 3 , and A = 0.18 eV, resulting in 
the spin-flip time r s = l/(A/3) = 7.7 x 10~ 13 s. Note 
that the spin flip length for permalloy is l s ~ 5 nm (44| 
corresponding to a spin flip time of t s — 2.3 x 10~ 14 s. 
A spin current density of J s = 0.2 then corresponds to 
a charge current density (with polarization P = 1) of 
I/A = 7 x 10 12 A/m 2 . Typical current densities realized 
in experiments are of this order of magnitude [33l . l34j . 

The numerical results for the DW configuration at 
times t = 80to and t = llOto for J s = 0.1 are shown 
in Fig. Q] As can be seen, the initial DW is moved to 
the right due to the presence of the spin current. Im- 
portantly, the second-order spin torque leads to a con- 
siderable deformation of the initial shape of the DW, 
where the magnetization rotates into the x — y plane 
around the center of the DW. This region grows with 
time, but the DW center still moves along the direction 
of the current. For comparison, we also depict the re- 
sults when only the first-order spin-torque is taken into 
account. This shows that the smearing of the DW profile 
is caused by the second-order torque. The strong DW 
deformation is the reason why a simple Walker Ansatz 
[|, H, H IH does not work. In the Walker Ansatz, 
one takes an adiabatic DW shape whose center is mov- 
ing with some velocity Vow(t), a spatially independent 
out-of-plane angle 4>{t), and a trial function for the po- 
lar angle 0([x— J d,TV£>w( T )]/w(t)) with time-dependent 
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FIG. 1: (Color online) DW configuration at times t = 80to 
[black solid line with filled circles] and t — llOto [red solid 
line] for J s =0.1 starting from the initial DW configuration 
(|34[) [dotted line] . For comparison, the result of keeping only 
Ti is also shown (for t = llOto, dashed line). 
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FIG. 2: (Color online) Time dependence of DW center for 
the same parameters as in Fig.[T] for J 3 = 0.05(»), 0.1(x) and 
0.2(D). Shown are the results with the spin torque included 
up to first-order (dashed lines) and up to second order (full 
lines). For J s = 0.05, both coincide in the shown regime of 
time (black lines and filled circles). 



DW width w(t). Inserting such an Ansatz into the full 
LLG equations yields four coupled differential equations 
for three variables (<j>(i),VDw(t) and w(t)). Unfortu- 
nately, there is thus no consistent solution under a Walker 
Ansatz. 

Despite the deformation of the DW, the DW center 
Xdw (defined as the zero of the n 2 -component) dis- 
plays a uniform steady-state translational motion. This 
is shown in Fig. [21 where Xx>\y(t) is shown for different 
J s . For the smallest value of the spin current density, 
J s = 0.05, the first- and second-order results for Xowit) 
coincide in the shown regime of time. For larger spin cur- 
rents, however, the DW center initially moves mainly due 
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0.25 




FIG. 3: (Color online) DW velocity vs. spin current density 
J s for same parameters as in Fig. [1] The black solid line gives 
the full numerical result, red empty squares the one when 
keeping only the first-order spin torque. The red dashed line 
gives the analytical result of Ref. [13| (for the full first-order 
spin torque), the blue full line with filled circles gives the 
result of Ref. [j| where only the spin torque is kept. 



to the first-order torque, and only after a transient time 
the DW is deformed. The motion of the DW center is 
then slower but still uniform, and a constant DW velocity 
can be extracted. At long times Vnw(t) = dXjjw(t)/dt 
approaches the constant value Vdw, shown in Fig. [3] as a 
function of J s . We also show the DW velocity following 
from only the first-order torque, where our numerical re- 
sults reproduce the known result Vow — flJ s /o> [H[. As 
one can see, for small J s , the first- and second-order re- 
sults merge. However, for growing J s , while the adiabatic 
(first-order) calculation gives a steadily growing DW ve- 
locity with growing J s , the second-order (non-adiabatic) 
calculation indicates a non-monotonic dependence. This 
is a central result of our work, as it clearly shows that a 
larger spin current does not necessarily increase the DW 
velocity. Instead, it may lead to a DW deformation and a 
slower DW steady motion. This finding seems in qualita- 
tive agreement with experimental observations where the 
measured DW velocities are typically smaller than those 
predicted by the so far available first-order calculations 
|33j ]. We note that deviations from the first-order result 
become noticeable around J s ~ 0.05 which corresponds 
to a charge current density oil /A = 1.75 x 10 12 A/m 2 for 
our parameters (with P = 1). This value should be com- 
pared, for instance, with the applied current densities of 
I/A = 1.0 xlO 12 A/m 2 inRef.iaand of I/A = 2.2 xlO 12 
A/m 2 in Ref. [34], indicating that the contribution of the 
second-order spin torque can be significant already for 
intermediate realistic current densities. 

In Fig. [3J we also show the result of Ref. 0] where 
only the spin torque JT]) was included, leading to Vow — 

\J J'j — Js cr w ith the critical threshold spin current 

Js,cr = K^w. The DW would then be pinned by the 
transverse anisotropy for small J s . However, as has been 



discussed in Ref. [HI , the full first-order spin torque lifts 
this pinning. 

Finally, we note that our formalism allows to numer- 
ically solve the three coupled equations of motion for 
n(x, t), s(x, t) and 3(x,t) directly. In particular, Eqs. (J3)), 
(fTTj) and (H2]), together with Eqs. dTHJ) and dTTJ), with the 
initial configurations given in Eq. (|34[) for n(x, 0) and 
s(x,0) = — x s An(a;,0) and J(x, 0) = J s n(x, 0), consti- 
tute nine coupled first-order ordinary differential equa- 
tions which can be integrated by a simple fourth-order 
Rungc-Kutta scheme. Since this scheme does not yield 
a simple closed formula for the effective spin torque T, 
however, we refrain from this approach here. 



VI. CONCLUDING REMARKS 

In this paper, we have studied the spin torque 
due to an applied spin-polarized current entering the 
Landau-Lifshitz-Gilbcrt equation for the magnetization 
— M s n(x, t) of a ferromagnetic metallic wire. Spin relax- 
ation of the itinerant electrons has been included within a 
Bloch-Redfield formalism adapted to this problem. The 
equation of motion of the itinerant electrons can be it- 
cratively solved within a gradient expansion scheme for 
current-induced spin torque T entering the LLG equa- 
tion. The first order in this expansion contains only 
first-order derivative terms (in particular, c^n), and re- 
covers the known 'adiabatic' result |T]) and the 'non- 
adiabatic' correction ((2]). We prefer to call both terms 
'adiabatic' here, since they are both of first order in the 
gradient expansion. The correction term discovered in 
Ref. [IH is in fact due to spin relaxation [IH, [H, 
We have then explicitly derived all second-order contri- 
butions to the spin torque term. As they involve spatial 
and time derivatives of the magnetization, they gener- 
ally induce a deformation of the magnetization profile as 
compared to the simple adiabatic shape. This has been 
illustrated for the example of a tail-to-tail transverse do- 
main wall. Despite the fact that the DW gets distorted, 
it still moves uniformly with constant terminal velocity. 
We have shown that this velocity depends in a non-trivial 
way on the applied spin current density J s . For small J s , 
it essentially coincides with the known first-order result, 
which illustrates that the 'adiabatic' pinning of the DW 
is overcome by the /3-term ([5]) in the spin torque. For 
larger spin currents, the DW velocity does not grow fur- 
ther with growing J s but instead diminishes again. Most 
likely, upon inclusion of even higher orders, this indi- 
cates a saturation of the DW velocity for large J s . Im- 
portantly, a simple analytical expression for the leading 
non-adiabatic corrections to the known current-induced 
spin torque T a d-\-TzL has been provided in our work, see 
Eq. (f3"2"| . We hope that this result will also find its sig- 
nificance in more elaborate micro- magnetic simulations. 
Finally, we remark that also the problem of pinning Q 
of a DW by structural defects could be included in our 
formalism. The generalization to incorporate thermal cf- 
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fects [46|, |47[ via a stochastic LLG equation appears pos- 
sible as well. 
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with fluctuating magnetic fields B lljq (x,t) produced by 
the harmonic bath. A key quantity is the (free) bath cor- 
relation function B aa (x, t)B VM i (x', t — r), or its Fourier 
transform, the spectral density [32|]. We will assume 
the spectral density to be real; neglecting its imagi- 
nary part only introduces small Lamb frequency shifts 
I28l. [29l. [30. [31I of no interest here. Thus we define 



APPENDIX A: RELAXATION KERNELS 

In this Appendix, we sketch the derivation of the re- 
laxation kernels entering the equations of motion for 
the spin density and spin current density of the itiner- 
ant electrons. We use the abbreviations K s (a;) = s(x) 
and Kj(i') = 3(x), and start from the master equation 
for the reduced density operator under the Markov ap- 
proximation for our time-dependent system Hamiltonian 
Hs(t) = Hq + H ex (t) . This master equation has been de- 
rived in detail in Ref. [30j | , and can be taken over directly. 
The derivation assumes that the system is linearly cou- 
pled to a collection of environmental harmonic oscillators 
(with 'bath' Hamiltonian Hb and bath temperature T) 
via a bilinear system-bath coupling Hs b ■ Specifically, we 
take 



H SB 



E 



q—x,y,z;X=s,J 



dxBx, q (x,t)K x , q {x), (Al) 



kqq',ij,u(x,x' , w) = - / dr cos(wt) 



(A2) 



x B ll , q (x,t)B„, q ,(x',t- t), 



which is independent of t. For our purpose, it is suffi- 
cient to assume a spatially independent isotropic spec- 
tral density which yields an x-independent correlator 
k q .^{w) = k qq i^ v (x,x' ,u)5 qq >. The coupling H SB is now 
assumed weak enough to justify a lowest-order perturba- 
tive scheme. Moreover, the bath autocorrelation time 
t c should be sufficiently short to justify the standard 
Markov (no memory) approximation. Our starting point 
is thus the Markovian master equation for the reduced 
density matrix p(t) of the system j30j , 



dp(t) 
dt 



-i[H s {t),p{t)]-- 



dT[H SB (t), [e-' H sW-H SB {t - r)e iff sW^p(t) - p s (t)]}, 



(A3) 



where the time-dependent Boltzmann density operator is 
given by 



equation of motion 



Ps(t) = 



e -H s (t)/k B T 
tTe -H s (t)/k B T- 



(A4) 



The system thus relaxes to the instantaneous Boltzmann 
distribution (|A4[) which contains the full time-dependent 
magnetic field generated by n(x,t) in H ex (t). 

From this result, we can then derive the Bloch-Rcdfield 
equations, i.e., the equations of motion for the quantum 
statistical expectation values (s(x,t)) and (J(x,t)). (For 
notational simplicity, we omit the brackets and the over- 
line (•) in most of the paper, but we always mean the 
quantum statistical expectation value. In this Appendix, 
we shall keep them for clarity.) Let us now consider the 
r-th component of (K\(x, t)) (A ~ s, J), which obeys the 



■(K x , r (x,t))=Y^ 



dpaa' (t) 



(a'\K x Jx,t)\a), (A5) 



CLOi' 

since there is no explicit time dependence present in the 
operator K>. We have employed a suitable complete 
set of cigenstatcs, Hs(t)\a(t)) = a(t)\a(t)), where for 
simplicity a denotes both the quantum numbers and 
the eigencnergy, and X]al a ){ a l = 1 ( we keep the t- 
dependence of a implicit from now on) . By inserting Eq. 
(|A3|) into Eq. (|A5[) , we obtain first a part simply yielding 
the coherent terms in Eqs. and (TT2"|) . Here we focus 
on the relaxation part and discuss the term ~ p(t) in 
Eq. (|A3[) . The other term involving ps(t) follows from a 
simple substitution. Therefore, we have to calculate 
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-(K x , r (x,t)) rel 



1 r°o 

/ dT(a\[HsB(t),[e- iH °> 



tyH SB (t-Ty H s(t)-, P {t))]\a>)(a'\K x , r (x)\a). (A6) 



r 



Expanding the double commutator yields four terms, of 
which one will be discussed in more detail. Inserting 
unities, we find such a term as 



dT(a\H S B(t)\f3){p>\H SB (t - T )\a') 



xe 



i(a' — f3')r 



{f3\p\P'){a'\KxAx)\a)- 



Using the bath correlation functions introduced above, 
this can be written as 

HZ! / dx'dx"{j3\K VA {x")\a) 

x{a\K x , r (x)K tJi]q (x')p\P}k q ^ u ((3 - a). 
Collecting all four terms then yields the relaxation part 
d 



(K X , r (x)) 



t.rel 



dt 

x(a\[[K x ^),K^ q (x%p]\l3)k q .^(f3~a), (A7) 

where we exploit that k q ^ v (uj) is even in to. 

Applying Eq. (|A7j) to the case K S)T=Z = s z is straight- 
forward upon using the Kac-Moody algebra and cyclic 
invariance of the trace. In the g-summation, the contri- 
bution for q = x involves the matrix elements of the oper- 
ator s x . States a and (3 joined by s x thus have (3 — a = A, 



the Larmor frequency for this problem. A similar reason- 
ing applies to q = y, while q = z does not contribute at 
all. Collecting all pieces yields for the relaxation part 



1 



with 
1 



-r;{sz(x,t)} rel = (s z (x,t)), 

dt t sz 



k x . ss {A) + k ViSS (A) + v 2 [k xJJ (A) + k yJJ (A)}, 



where we use k q ,^{uj) — k q , v ^{uj). For r = x,y, we 
proceed in the same way and find 



1 



with 
1 

T s,x/y 



^-(S x/ y(X,t)) rel = -—^{S X /y{X,t)) 



k y/x ^ s (A) + k z , ss {0) + v 2 [k v/xJ j(A) + k z jj(0)}. 



Assuming isotropy, t s _ x « t s>v « r s z w r s , and adding 
the part from the stationary distribution (|A4[) . we finally 
obtain the relaxation term (| 16[) entering Eq. (fTTj) . In the 
same manner, Tj follows from such a calculation. 
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